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AHHOTan,HH 

M.IO.BorojiK)5cKHH h pp. KoppeKii,Ha HejiHHeiiHOCTH SHepreTHnecKoii niKajiBi sjieKTpoMarHHTHoro KajiopHMeTpa 
no flByxcJaoTOHHbiM pacnaflaM Tr^-Meaona.: HpenpHHT HOBB 2010-15. - HpoTEHHO, 2010. - 6 c, 4 pnc, 6h5- 
jiHorp.: 4. 

B pa5oTe npeflCTaBjien MeTOfl BbiHHCjieHHa KoppeKn,HH HejiHHeitHOCTH OTKjiHKa ajieKTpoMarHHTHoro Kajiopn- 
MBTpa, ocHOBaHHbiH Ha MHHHMH3aii,HH OTKjiOHeHHH H3MepeHHolt Maccbi HeHTpajibHoro MesoHa, pacnaflaMmeroca 
B KOHeHHOM CHeTB Ha (J'OTOHbi, B saBHCHMOCTH SHeprHH HOCjieflHHx. MeTOfl 5bijT pa3pa6oTaH H npHMenen fljia 
3JieKTpoMarHHTHoro KajiopHMexpa LGD2 b SKcnepHMeHTe FHnepoH-M na ycKopHTejie YTO FHIl HOB3. Haii- 
flCHHaa KoppeKii,Ha nosBOjiHjra cymecTBeHHO yMeHbniHTb BapHaLi,HH peKOHCxpyHpoBaKHbix Mace tt" h 77 MesoHOB b 

SaBHCHMOCTH OT HX MHHHMajIbHOH SHeprHH. 

Abstract 

M.Yu.Bogolyubsky et al. Correction of the energy scale nonlinearity in electromagnetic calorimeters with the tt" 
two-photon decays: IHEP Preprint 2010-15. - Protvino, 2010. - p. 6, figs. 4, refs.: 4. 

The method to calculate the non-linearity correction of the electromagnetic calorimeter response, based on 
minimisation of the deviation of the measured neutral meson mass on the energies of it decay photons, is described 
in this paper. This method was developed for the electromagnetic calorimeter LGD2 in the Hyperon-M experiment 
at U70 accelerator of IHEP. The found correction allowed to reduce significantly variations of the reconstructed 
7r° and 77 masses on the minimal energy of the mesons. 
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Introduction 



Photons and electrons due to interaction with a medium of the cell-type electromagnetic calorimeter 
produce electromagnetic showers which spreads over several calorimeter cells called a shower cluster, i.e. 
the group of affected cells with common edges. The read-out electronics for such kind of calorimeters 
reads the signal amplitudes from calorimeter cells. These amplitudes are used to estimate the real 
energy deposition of electromagnetic shower in the calorimeter cells by using the independent on energy 
calibration coefficients. The sum of the deposited energies in the cluster cells defines the energy of the 
incident photon or electron. This direct energy estimation of electromagnetic showers might be satisfactory 
in the energy range used for the calorimeter calibration but could lead to energy shifts at different 
energies which results in the calorimeter response nonlinearities caused by the physical processes, read- 
out electronics and shower reconstruction program. 

The longitudinal electromagnetic shower profile (electromagnetic cascade in the calorimeter radiators) 
[T] allows to determine the shower energy deposition in the calorimeter radiators for the case of its finite 
longitudinal thickness. The position of the energy maximum moves further into the calorimeter with 
the logarithm of the photon energy, that increases the shower energy leakage out of the calorimeter. 
Another phenomenon of the measured shower energy loss is related to the finite attenuation length for 
Cherenkov or scintillation light in the calorimeter cells. The average light path from a radiation point to 
a photo-detector depends on the energy of the incident photon and reveals itself also as the nonlinear 
dependence with energy of the light pulse produced by shower. The shower energy leakage is possible in 
the transversal directions as well, for instance, due to energy loss in gaps between calorimeter cells. 

Chosen calorimeter design could bring the nonlinearity effects as well. For instance, the used photo- 
detectors could have a nonlinear scale. The read-out electronics (including the analog to digit converters, 
ADC) could be too noisy, and the noise has to be suppressed by applying the relevant threshold on 
recorded amplitudes in the calorimeter cells. This threshold leads sometimes to a significant distortion 
of measured amplitudes of the incident photons at low energies. The enumeration could be continued. 
But it is important to note that all these effects are unlikely possible to take into account with a high 
accuracy using Monte Carlo simulations only. Anyway this is sufficiently difficult. 

The typical task solving by electromagnetic calorimeters in high energy physics experiments is the 
mass spectra measurement of neutral mesons decaying into photons, for instance, tt" 77, 77 77, 
ui — >■ 7r°7 and so on. The calorimeter energy scale nonlinearity have an impact on dependence of the 
measured neutral meson masses on their energies which leads, in turn, to systematic uncertainties in the 
meson spectra measurement. Therefore the correction of the calorimeter non-linearity response is relevant 
in the case. 

At the same time the possibility of solving this problem directly is no means always the case, i.e. the 
experimental study of the calorimeter response to photons or electrons at different energies cannot be 
carried out, for example, at collider experiments or for other reasons. However the correction factor of 
energy scale of electromagnetic calorimeters could be found as the result of inverse problem solution, i.e. 
by using experimentally measured mass dependence of neutral mesons on the energy of decay photons. 
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In the present paper the mathematically strict algorithm of nonlinearity correction of the calorimeter 
energy scale based on the minimum squared deviation method is proposed. This algorithm has been 
developed and applied for the data processing from the electromagnetic calorimeter LGD2 of the experiment 
Hyperon-M at the U-70 accelerator of IHEP, Protvino [2]. The two photon decays of neutral pions recorded 
in the experiment have been used for the energy correction procedure. The performed correction allows to 
reduce significantly the nonlinearity of the LGD2 energy scale and to decrease systematic uncertainties in 
particle mass measurement in several times. It opens up the possibility to obtain the interesting physics 
results as well. 

It is worth to note also that the events of two photon decays of neutral mesons are used for a calibration 
purpose of the relevant electromagnetic calorimeters in several experiments. And thus the described below 
procedure of the energy scale correction could be interesting for the data treatment in these experiments 
as well. 

1. Experiment Hyperon-M 

It is appropriate at first to give a short description of the Hyperon-M setup before discussing the 
electromagnetic calorimeter LGD2 energy scale in the experiment. The layout of experiment is presented 
in FigU] The setup comprises the beam telescope of scintillation counters Si, S2, S4, Cherenkov counters 
Ci_3, nuclear target T, scintillation anti counter Sa and electromagnetic Cherenkov lead glass calorimeter 
LGD2 located at a distance of 3.7 m after the target. The measurements were carried out on the 7 GeV/c 
beam of positive particles with intensity of 10® particles per burst on different nuclear targets, including 
the Be target. The requirement of a beam particle signal from the beam telescope and the absence of a 



LGD2 




Phc. 1. The Hyperon-M experimental setup layout: Si, S2, S4 — beam scintillation counters, Ci_3 — Cherenkov 
counters, T — nuclear target, Sa — trigger scintillation anti counter Sa, PCi — proportional chambers, 
LGD2 — Cherenkov electromagnetic calorimeter with lead glass radiators. 

signal from anti counter Sa generates the trigger signal: 

Tr = Si - 82- Sa - Sa- 

This trigger allows to select effectively the inclusive production of neutral mesons decaying into 
photons within the LGD2 calorimeter solid angle: 

TT+{K+,p) + A^^ M" + X, Af° ^ n7. (1) 

A typical value of the trigger selectivity reached the value of ^ 1-3 • 10^^ depending on the type and 
thickness of the irradiated target and the beam intensity. More detailed description of the Hyperon-M 
setup, electronics, trigger and data acquisition system can be found elsewhere [3]. 

The LGD2 calibration was performed on the physics two-photon events collected on the Be target. The 
sample of calibration events comprises of 2 millions events ([1]) with the reconstructed photon multiplicity 
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n = 2 and the photon pair energy E^^ > 1.5 GeV. Determination of the cahbration coefficients was 
performed by means of the iterative corrections of the Tr'^-peak position with a smooth background in 
each calorimeter cell on the subset of two photon events where one of two photons hits this cell, details 
see in [4]. We note here only that the effective mass of photon pair was evaluated with the formula: 

m2-y = v/2£ie2(l - cos 6*12), (2) 

where Si is the measured energy of the i-th photon and 9i2 is the opening angle of photon pair in 
the laboratory frame. The effective mass spectrum of photon pairs in reaction ([1]) after 15 iterations is 
illustrated by Fig. ^ The obtained mass resolution for the 7r°-meson is equal to 11.4 MeV. 



2. The calorimeter energy scale correction procedure 

Let's define the nonlinear correction to the calorimeter LGD2 energy scale Ae as the difference between 
the "true" photon energy e and its measured value e: 

Ae^e-e. (3) 

This correction can be expanded in a power series over some variable x depending on the photon energy 

i—k 

Ae = ^ai-x' (4) 

i=0 

taking into account that the correction Ae should be comparatively small with respect to the measured 
photon energy. To avoid the computational precision limitations at large energy values related to the 
factorisation order k in expression (0]), it is natural to take for the x variable the logarithm of measured 
photon energy: 

X = x{e) = ln(e/£o), (5) 
where £0 = 1 MeV. As a consequence the corrected photon energy e can be written as: 



i—k 

e{e)^e + Ae^e{l + J2—^l, (6) 
i=0 ^ 



where it is natural to assume that the parameters ai/e are sufficiently small due to a small nonlinearity 
of the calorimeter energy scale. The expression for the effective mass of a photon pair ^ can be rewritten 
then in terms of the corrected energies of photons as follows: 



TO27 = \/2£i£2(l - cos 6*12) = y^l^-Ci2, (7) 



where ii = £(£i) are linear functions ^ of small parameters ai/e and C12 = vT— cos5i2 is the geometrical 
factor which is actually independent on these parameters. 

The parameters in equation ^ can be determined by minimisation of the deviation of effective 
mass of the photon pair in representation ([7]) from the PDG 7r°-meson mass on the sample of 7r° events 
used in the discussed procedure and shown for our case in Figl^] (left) as hatched area. In other words, 
the parameters Oii can be determined by means of the functional minimisation 



N ,~ ^2 

X = 1^ ^2(„,^^) ' (8) 



=1 



where N is the number of two-photon events in the indicated Tr'^-peak region in FigUl m27 is the effective 
mass of a photon pair in the representation ([7]), nij^o is the PDG value of the Tr^-meson mass and (7(77127) 
is the expected uncertainty of the effective pair mass as defined in expression 

The uncertainties on the invariant mass of photon pair include the photon energy uncertainty and the 
uncertainty of the photon pair opening angle, see ([2]). The opening angle error is defined by the Hyperon- 
M setup geometry and the reconstruction program of LGD2 calorimeter. This error is sufficiently small 
in our case, and we will neglect it below. 
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The relative uncertainty of the photon energy measurement in electromagnetic calorimeter is defined 
according to the formula: 

cTe/e = a/v^® 6 ® c/e, 

where parameters a, b and c are defined by the calorimeter design, see for example [1]. The last summand 
contribution in the energy resolution of LGD2 calorimeter is small and we will ignore it below. Thus the 
expected mass resolution for photon pairs can be express using the error propagation techniques as: 

a\m2y) = A{cl^{ei+e2) + B), (9) 

where the energies of photons are measured in GeV, and A and B are the empirical parameters equal to 
2.5 • 10"^ GeV and 1.4 • 10~^ GeV respectively for the LGD2 spectrometer. 
The necessary conditions for functional minimisation 

dx'^/da, ^ 

with the accuracy up to the second order smallness a^aj/e^ result in the system of linear equations 
relatively to the parameters aj: 

k N 2 

53 . ^2(^ \ (^1^2 + e2x\){eix{ + £23^1) = (10) 

N 

= V(^^ - Ci2)^^^{eix\ + £2X1), 



where xi = x{ei), / = 1,2, see equation ©. The iteration procedure based on equations ^ and ([TU]), 




MeV e.^, GeV 

Phc. 2. Left: the eflFective mass spectrum of two-photon events with the sum energy of photons larger than 1.5 
GeV. Right: the energy scale correction function for LGD2 calorimeter defined on the 7r° event sample 
and shown as hatched area in picture on left. 

allows one to find out the functional minimum ([5]) with a reasonably good accuracy after 2 — 3 iterations. 

The first 9 terms of the series (H]), i.e up to the order of fc = 8, were taken into account for the 
Hyperon-M data treatment. The next values for the correction coefficients have been obtained after the 
first iteration: ao-s/GeV = 0.00399, -0.0505, -0.0392, -0.0209, -0.00537, 0.0165, 0.0104, 
-0.00313, -0.00235. The coefficients of the second iteration were found to be about three times less 
compared with the first iteration values. The energy correction function for the LGD2 calorimeter 

i—k 

e™^7e = £{£)/e = l + ^ ^ (11) 
i=0 ^ 
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is presented in Figl^] in the right plot. As one can see from the figure the relative correction doesn't 
exceed 10% level virtually in the whole energy range of photons and this is in a good agreement with 
our initial assumption concerning the smallness of an expected energy scale nonlinearity of the LGD2 
calorimeter. This is an important statement because it is used in the ground of the method. For the sake 
of completeness it would be useful also to present the values of ® before and after the correction: in 
our case the value of per degree of freedom before the correction and after it are equal to 1.073 and 
1.044 respectively for approximately 10^ degrees of freedom. 

3. Results and discussion 

Performance of the above discussed procedure is illustrated in Fig[21 where the scatter plots of the 
effective two photon mass versus the logarithm of the energy of each photon in a pair (two points per event) 
for the reconstructed two-photon events ([T]) is shown for _Be-target before the energy scale correction on 
the left panel and after it on the right panel of the figure. A clear correlation of the two-photon mass 
and the photon energies for events in the Tr^-meson region is seen on the left picture and it is completely 
absent on the right one. The numerical values of the correlation coefficients for events without the energy 
scale correction and with it are equal to 0.13 and 0.05 respectively. 




ln(e/eQ) ln(e/eQ) 

Phc. 3. Distributions of the two-photon effective mass versus the logarithm of photon energy of each photon in 
pair (two points per event) for the reconstructed two-plioton events on Be-target. The concentration of 
events at lower area corresponds to the detection of Tr^-mesons, upper one - to the jy-mesons. The event 
distributions before the energy scale correction are shown on left and after the correction - on right. 

Another illustration of the nonlinearity correction method is represented by Fig. |4l The dependence 
of the measured mass of the tt^- and r^-mesons on the minimal photon pair energy (ffsum = ei + £2) is 
shown before applying the nonlinearity correction and after it. These plots demonstrate as well that the 
systematic deviation of the neutral pion mass from the PDG value in dependence on the photon pair 
energy decreases from 1.17% to 0.19%, i.e. in 6 times, and the same deviation for ?7-meson decreases 
from 2.98% to 0.23%, i.e. in 13 times, and this is demonstration of the high performance of the proposed 
method as a whole. 

Conclusion 

This paper describes the procedure of the energy scale correction for electromagnetic calorimeters. 
The procedure is based on the minimization of mass resolution for two-photon decays of neutral pion 
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Phc. 4. The mass dependence of 7r°-meson (on left) and jy-meson (on right) versus the minimal photon pair 
energy Ssum = ei + £2- The dependence before the LGD2 energy scale correction is shown by black 
colour and one after correction is shown by red colour. The dashed line shows the PDG values of these 
mesons. 



detected in the calorimeter. The linear parametrisation of the correction function as the power series in 
logarithm of the photon energy allows to provide a simple and effective energy scale correction in a very 
wide energy range. Possibility to use the physics statistics of the experiment for the energy correction 
procedure results in the high accuracy and sensitivity of the method. For instance, in the Hyperon-M 
experiment the reached mass scale nonlinearity for two-photon events is equal to 0.2%, that hardly can 
be obtained in calculations of similar corrections by Monte-Carlo methods due to restrictions peculiar to 
the transport code. Anyway the significant Monte-Carlo difficulties appear in calculations at the accuracy 
level of 10~3. 

Eventually, it is significant that the described procedure could be applied practically for any hodoscopic 
electromagnetic calorimeter if the physics statistics of experiment possesses the needed amount of two- 
photon or, let's say, three-photon decays of known mesons because this procedure could be easy generalised 
for multi-photon decays as well. 

The authors appreciate N.A.Kuzmin for discussions and helpful comments. 
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